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O ' Abstract 

The determination of the densest packings of regular tetrahedra (one of the five Platonic solids) 
is attracting great attention as evidenced by the rapid pace at which packing records are being 
broken and the fascinating packing structures that have emerged. Here we provide the most general 

o 

^D \ analytical formulation to date to construct dense periodic packings of tetrahedra with four particles 

^^ ■ per fundamental cell. This analysis results in six-parameter family of dense tetrahedron packings 

o: 

that includes as special cases recently discovered "dimer" packings of tetrahedra, including the 
densest known packings with density (j) = |2§y = 0.856347. . .. This study strongly suggests that 
^ • the latter set of packings are the densest among all packings with a four-particle basis. Whether 

they are the densest packings of tetrahedra among all packings is an open question, but we offer 
remarks about this issue. Moreover, we describe a procedure that provides estimates of upper 
bounds on the maximal density of tetrahedron packings, which could aid in assessing the packing 
efficiency of candidate dense packings. 

PACS numbers: 61.50.Ah, 05.20.Jj 



I. INTRODUCTION 



Dense packings of nonoverlapping solid objects or particles are ubiquitous in synthetic 
and natural situations. Packing problems arise in technological contexts, such as the pack- 
aging industries, agriculture (e.g., grains in silos), and solid-rocket propellants, and underlie 
the structure of a multitude of biological systems (e.g., tissue structure, cell membranes, 
and phyllotaxis). Dense particle packings are intimately related to the structure of low- 
temperature states of condensed matter, such as liquids, glasses and crystals [l^^ . In the 
last decade, scientific attention has broadened from the study of dense packings of spheres 
(the simplest shape that does not tile Euclidean space) |5l-ll4| to dense packings of disor- 
dered 15Nl9| and ordered 20N22| nonspherical particles. In addition, the equilibrium phase 



behavior and transport properties of hard nonspherical particles have been of topics of great 
interest 



23 



26| . Nonsphericity introduces rotational degrees of freedom not present in sphere 



packings, and can dramatically alter the characteristics from those of sphere packings. 

A basic characteristic of a packing in rf- dimensional Euclidean space M.'^ is its density 0, 
defined to be the fraction of Mf^ that is covered by the particles. A problem that has been 
of great scientific interest for centuries is the determination of the densest arrangement (s) 
of particles that do not tile space and the associated maximal density (pmax- For generally 
shaped particles, finding the densest packings is notoriously difficult. This salient point is 
summarized well by Henry Cohn |27| who recently remarked, "For most grain shapes we 
cannot guess or even closely approximate the answer, let alone prove it, and it is difficult to 
develop even a qualitative understanding of the effects of grain shape on packing density." 
Until recently, very little was known about the densest packings of polyhedral particles. The 
difficulty in obtaining dense packings of polyhedra is related to their complex rotational 
degrees of freedom and to the non-smooth nature of their shapes 
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29|. 



Recently, we set out to attempt to determine the densest known packings of the Platonic 
and Archimedean solids l28l 



29| . It was shown that the central symmetry of the majority 



of the Platonic and Archimedean solids distinguish their dense packing arrangements from 
those of the non-centrally symmetric ones in a fundamental way. (A particle is centrally 
symmetric if it has a center C that bisects every chord through C connecting any two bound- 
ary points of the particle; i.e., the center is a point of inversion symmetry.) The tetrahedron 
is the only Platonic solid that lacks central symmetry, an attribute that geometrically frus- 



trates it to a greater degree than the majority of the remaining sohds in this set that do 



not tile space 29|. A number of organizing principles emerged in the form of conjectures for 
polyhedra as well as other nonspherical shapes. In the case of polyhedra, the following three 
are particularly relevant and will be applied later in the paper to remark on the optimality 
of tetrahedron packings: 

• Conjecture 1: The densest packings of the centrally symmetric Platonic and 
Archimedean solids are given by their corresponding optimal Bravais lattice packings. 

• Conjecture 2: The densest packing of any convex, congruent polyhedron without 
central symmetry generally is not a Bravais lattice packing, i.e., the set of such poly- 
hedra whose optimal packing is not a Bravais lattice is overwhelmingly larger than the 
set whose optimal packing is a Bravais lattice. 

• Conjecture 3: The densest packings of congruent, centrally symmetric particles that 
do not possesses three equivalent principle axes (e.g., ellipsoids) generally cannot be 
Bravais lattices. 

Conjecture 1 is the analog of Kepler's sphere conjecture for the centrally symmetric 
Platonic and Archimedean solids. In this sense, such solids behave similarly to spheres 
in that their densest packings are lattice arrangements and (except for the cube and the 
truncated octahedron) are geometrically frustrated like spheres. Conjecture 2 has been 



shown by Conway and Torquato 30|] to be true for both the tetrahedron and Archimedean 



non-centrally symmetric truncated tetrahedron, the latter of which can be arranged in a 
"uniform" non-Bravais lattice packing with density at least as high as || = 0.958333 .... (A 
uniform packing has a symmetry operation, in this case the point inversion symmetry, that 
takes any particle into another.) Although Conjecture 3 is explicitly stated in its current 
form here for the first time, it was strongly implied in Ref. 29 |. 



It was Conway and Torquato's investigation [30| that has spurred the flurry of activity 
in the last several years to find the densest packings of tetrahedra. There have been many 
twists and unexpected turns since 2006 that have led to the dense packings of tetrahedra 
that we report here. Therefore, to place our present results in their proper context, it is 
instructive to review briefly the developments since 2006. 



TABLE I: A brief summary of the dense non-lattice packings of tetrahedra. The name of the 
packing is given along with the year that it was discovered. Here (p is the packing density and N 
is the number of tetrahedra per fundamental cell. 



Packing 




Year 





N 


Uniform I [30] 




2006 


1 = 0.666666 . . . 


2 


Welsh [30j 




2006 


1| = 0.708333... 


34 


Icosahedral [30] 




2006 


0.716559... 


20 


Uniform II [35] 




2009 


i^^fgV^o -0.719488.. 


2 


Wagon Wheels [33] 




2008 


0.778615... 


18 


Improved Wagon Wheels 


[28] 


2009 


0.782021 . . . 


72 


Disordered Wagon Wheels 


[29] 


2009 


0.822637... 


314 


Ring Stacks [34] 




2009 


0.8503... 


82 


Uniform III [35] 




2009 


^ = 0.854700 . . . 


4 


Dimer-Uniform I [36] 




2009 


iflfl = 0.855506 . . . 


4 


Dimer-Uniform II [37] 




2010 


^°°? -0.856347... 


4 



First, we note that the densest Bravais-lattice packing of tetrahedra (requiring one tetra- 
hedron per fundamental cell such that each tetrahedron in the packing has the same orien- 
tation as the others) has a packing fraction = ;|| = 0.367 . . . and each tetrahedron touches 
14 others 3l|] • Conway and Torquato [30] showed that the densest packings of tetrahedra 
cannot be Bravais lattices by analytically constructing several such packings with densi- 
ties that are substantially larger than ||. (A non-Bravais lattice packing contains multiple 
particles, with generally different orientations, per fundamental cell, which is periodically 
replicated in M^.) One such packing is a "uniform" packing with density = 1 and two 
particles per fundamental cell. The so-called "Welsh" packing has a density (j) = 0.708333 . . . 
and 34 particles per fundamental cell. Yet another non-Bravais lattice packing with den- 
sity = 0.716559 ... is based on the filling of "imaginary" icosahedra with the densest 
arrangement of 20 tetrahedra and then arranging the imaginary icosahedra in their densest 
lattice packing configuration. The densities of both the Welsh and Icosahedral packings 
can be further improved b y c ertain particle displacements 30|. Using imperfect "tetrahe- 
dral" dice, Chaikin et al. J32| experimentally generated jammed disordered packings with 



~ 0.75. Employing physical models and a computer algebra system, Chen [33| discovered 
a remarkably dense periodic arrangement of tetrahedra with cj) = 0.7786 . . ., which exceeds 
the density {(pmax = tt/vIS = 0.7404 . . .) of the densest sphe re p acking by an appreciable 
amount. We have called this the "wagon- wheels" packing 28|, |29j. 

Torquato and Jiao [28| devised and applied an optimization scheme, called the adaptive- 
shrinking-cell (ASC) method, that used an initial configuration based on the wagon-wheels 
packing to yield a non-Bravais lattice packing consisting of 72 tetrahedra per fundamental 
cell with a density = 0.782 . . . [28]. Using 314 particles per fundamental cell and start- 
ing from an "equilibrated" low-density liquid configuration, the same authors were able to 



improve the density to = 0.823 . . . 29||. This packing arrangement interestingly lacks long- 



range order. Haji-Akbari et al. 3^ numerically constructed a periodic packing of tetrahedra 
made of parallel stacks of "rings" around "pentagonal" dipyramids consisting of 82 particles 
per fundamental cell and a density = 0.8503 . . .. More recently, Kallus et al. 35|| found 
a remarkably simple uniform packing of tetrahedra with high symmetry consisting of only 
four particles per fundamental cell with density cj) = j^ = 0.854700 .... We subsequently 
presented an analytical formulation to construct dense uniform dimer packings of tetrahedra 
and employed it to obtain a three-parameter family of packings. (A dimer is composed of 
a pair of regular tetrahedra that exactly share a common face. A uniform dimer packing of 
tetrahedra takes any dimer via a point-inversion symmetry operation into another.) Making 
an assumption about one of these parameters resulted in a two-parameter fami 



those with density as high as 



12250 
14319 



y, including 
0.855506... [36|. Chen et al. 371 recognized 



that we made such an assumption and employed a similar formalism to obtain a three- 
parameter family of tetrahedron packings, including the densest known dimer packings of 
tetrahedra with a density = ^^ = 0.856347. . .. Table [T] summarizes some of the packing 
characteristics of the non-Bravais lattice packings of tetrahedra. 

In the following section, we provide the details of our more general formulation and 
construct a six-parameter family of dense tetrahedron packings. We will show that our 
brmalism includes as special cases all of the recently discovered four-particle basis packings 



35 



37|. Our analysis strongly suggests that the optima among this set of packings provide 
the densest arrangements among all packings with a four-particle basis. In Sec. IIIIl we 
describe a procedure that provides estimates of upper bounds on the maximal density of 
tetrahedron packings, which could aid in assessing the packing efficiency of candidate dense 



packings. In Sec. llVt we make concluding remarks, including comments on the optimality 
of the densest known dimer packings of tetrahedra. 

II. ANALYTICAL CONSTRUCTIONS OF DENSE PACKINGS OF TETRAHE- 
DRA 



Inspired by the work of Kallus et al. 35|, we have applied the adaptive- shrinking- cell 



(ASC) optimization scheme to examine comprehensively packings with a considerably small 
number of particles per fundamental cell (from 2 to 32) than we have used in the past 



28 



|29|. The ASC scheme employs both a sequential search of the configurational space 
of the particles and the space of lattices via an adaptive fundamental cell that deforms 
and shrinks on average to obtain dense packings. A dense packing with 8-particle basis 
that emerged from this numerical investigation suggested that it was composed of two very 
similar fundamental cells, each containing 4 particles. Using one of the 4-particle basis 
configurations, we were able to find packings with density = 0.8551034 . . . that exceeded 
the highest density packings with (p = ^ = 0.854708 . . . constructed by Kallus et al. Even 
though our packings possess a type of point inversion symmetry, they are not as symmetric 



as the densest packings reported in Ref. [35|, as we now explain. 

The four tetrahedra in the fundamental cell in our dense numerically generated packings 
formed two contacting "dimers" . A dimer is composed of a pair of regular tetrahedra with 
unit edge length that exactly share a common face. The compound object consisting of 
the two contacting dimers possesses point inversion symmetry, with the inversion center 
at the centroid of the contacting region on the faces. A Bravais lattice possesses point 
inversion symmetry about the lattice points and the centroids of the fundamental cells. By 
placing the symmetry center of the two-dimer compound on the centroids (or the lattice 
points), we construct packings that generally possess point inversion symmetry only about 
the symmetry centers of the two-dimer compound. Besides the centroids of the fundamental 
cell, all of the half-integer lattice points are also inversion symmetry centers of the packing. 
We call such structures dimer-uniform packings, since the inversion symmetry acts to take 
any dimer into another. Such packings should be distinguished from the more symmetric 
uniform (or transitive) packings of tetrahedra in which the symmetry operation acts to take 



any tetrahedron into another, such as the ones found in Refs. 30| and 35|| (see Table [T|). 



The latter have almost as much symmetry as a Bravais lattice, except that the centroids of 
the particles are not just characterized by simple translational symmetry. 

A. General Formalism: Six-Parameter Family of Dense Tetrahedron Packings 

We now describe our general analytical formulation to construct dense tetrahedron pack- 
ings by relaxing the symmetry conditions on the contacting dimers in detail. In particular, 
we orient the 3-fold rotational symmetry axis of one of the dimers in an arbitrary direction 
(say the z-direction of a Cartesian coordinate system), and then fix the origin of the lattice 
vectors at the centroid of this dimer. (The centroid is located at the center of the contacting 
faces of the two tetrahedra that comprise the dimer.) Then we place the second dimer in 
contact with the first one such that there is a center of inversion symmetry that takes one 
dimer to the other, which implies a face-to-face contact between the two dimers. 

The problem of determining the analytical constructions then amounts to determining 
12 equations for the 12 unknowns. Nine of the 12 unknowns arise from the three unknown 
lattice vectors, each of which contains three unknown components. The other 3 unknowns 
derive from the components of the centroid of the second dimer. 

In particular, we let the centroid of the dimer at the origin be denoted by Cq = (0, 0, 0), 
and the centroid of the other dimer be denoted by Ci = {'r]i,f]2,V3)- The vertices of the 
two dimers associated with Cq and Ci are given by v^ = (I, ^7^,0), v^ = {~2^ T7^^^)^ 

v^ = — vc + Ci, v|) = — V£) -|- Ci, v^ = — v^; + Ci, respectively. In addition, let the 
lattice vectors be Ai = (-ai, -0:2, -as), A2 = (/3i, -/32, -/^s), and A3 = (71,72,73)- The 
12 components of the four vectors Ci and Aj, {i = 1, 2, 3) are the variables that determine 
the packing. Note in the above general set-up, we assume no particular symmetry of the 
packings. 

In our packings, the dimers have 8 possible face-to-face contacts, 6 possible edge-to-edge 
contacts and 4 possible vertex-to-face contacts. Each contact provides a condition that 
either reduces the number of variables for the packing or constrains their feasible values. 

A face-to-face contact requires that the projection of the vector distance between the cen- 
troids of the two dimers on the contacting face normal is equal to the diameter of the insphere 
(the largest possible sphere that can be inscribed into a tetrahedron) of the tetrahedra. A 




FIG. 1: (color online). The geometrical set up for the general formulation of the analytical con- 
structions of dimer packings of tetrahedra. The dimer (two tetrahedra sharing a common face) 
centered at cq is shown in blue with the vertices labeled (note that v^; whose projection coincides 
with V£) is not shown). The dimer centered at ci is shown in red. The dimer centered at (cq + Ai) 
is shown in yellow. The dimer centered at (cq + A2) is shown in green. The dimer centered at 
(cq + A3) is shown in purple. Observe that in the perspective shown the dimers appear as if they 
were single tetrahedra. 

nonoverlapping condition associated with a face-to-face contact requires that the projection 
of the vector distance between the centroids of the two dimers on the contacting face normal 
is greater than the diameter of the insphere. The 8 possible face-to-face contacts are between 
the dimer pairs with the centroids at {ci,Co}, {ci,Co + Ai}, {ci,Co + A2}, {cq + A3,Ci}, 
{co + A3, ci - A2}, {co + A3, ci + Ai - A2}, {co, ci - Ai + A3}, and {cq, Ci + Ai - A2 - 2A3}. 
The possible contact between dimer pairs at {uj,Uj} requires 

2^6 



[u, - Uj) ■ riy > g 



(1) 



where rij^ is unit outward contacting face normal of the dimer at Uj. In the following 
discussion, we explicitly provide these equations. 

First, we note that Ci can be completely determined by considering the contacts that 



8 



the diiner associated with it makes with three other diniers centered at Cq, (cq + Ai) and 
(cq + A2) (see Fig. [1]). Note that Ai and A2 are two of the lattice vectors. The equations 
associated with these three contacts are respectively given by 

1 1 V2 1 ,,, 

:(?72 + a2) - -— ^(r/s + tta) = 77, (3) 



3x/3" ' 12^3'' ' 9 

By solving Eqs. ([2]) - @, we can obtain the components {rji,rj2,rj-i) of Ci, i.e., 

r/i = ^(6/3i- 2^3/32 -v^/33), (5) 



^^ 12 



-(-8a2 + 272^3 + 273/31 -2/32- 72/33), (6) 



^3 



-^ (4^3 - 4^2 + 72^3 - 2^3/31 + 2/32 + V2/33) , (7) 



3^2 
Moreover, we note that A3 can be completely determined by considering the contacts 

between the dimer centered at (cq + A3) and at Ci, (ci — A2) and (ci + Ai — A2), respectively. 

The equations associated with these three contacts are given by 

^(^1 - 71) - ^(^2 - 72) + ^(r^3 - 73) = \. (8) 

:(^2 + /32-72) + — ^(r?3 + /33-73) = 7T' (9) 



3v^" ' 12v^" '' 9^ 



1 1 /o 1 

--(^1 + «! - /3i - 7i) - 7-7^(^2 + a2 + 1^2- 72) + 777^(^3 + 03 + /33 - 73) = -. (10) 

By solving Eqs. ([8]) - f fTOj) . we can obtain the components (71,72,73) of A3, i.e., 

71 = ^(-6ai-2y3«2 + v^«3- 2^6/33), (11) 



^' 12 



^(-2y3ai- 10^2 + 3^203 + 8/32), (12) 



73 



-(-4v^ + V6ai + 3V2a2 - 2a^ + 2^6/31 + /^s 



(13) 



By achieving the aforementioned 6 face-to-face contacts, we can reduce the number of 
independent variables for the packing from 12 to 6, with the remaining variables being a^ 
and Pi {i = 1, 2, 3). In other words, all of the lattice vectors and the centroids of dimers are 
completely specified by the six parameters «« and f3i, which gives the six-parameter family 
of dense packings, i.e.. 



Cl 



^(6/3i-2v^/32- 76/33) 
ia2 + 2V2as + 2v^/3i - 2/32 - v^/^s) 
-^(4v^ - 4a2 + V2a3 - 2^3(3^ + 2/^2 + v^/^s) 



-1 T 



^< 



Ai = (-«! 



"2, 



tts), A2 = (/3i, -(32, -(33), 



(14) 



A, 



3^(-6«i - 2VSa2 + VQas - 2VQ/33) 
3^(-2V3ai - 10^2 + ^V2a3 + 8/32) 



i(-4v^ + Veai + 3V2a2 - 2^3 + 2v^/3i + (33) 
Note that these parameters can not be varied completely independently of each other, i.e., 
they are related by the additional nonoverlapping conditions which will be given in the 
ensuing discussions. 

There are two additional possible face-to-face contacts between dimer pairs centered at 
{cq, Cl — Ai + A3}, and {cq, Ci + Ai — A2 — 2A3}. The equations associated with these possible 
contacts are given by 



i(rh + a,+7i)-^ 



{r]2 + a2 + 72) 



V2 , ,1 

^r,3 + a3+73)>g, 



(15) 



6 



(771 - tti - /3i - 271) + — /^(^2 -a2 + (32- 272) 



^(r^3-a3 + /33-273)>^. 



(16) 



These conditions constrain the possible values of a^ and /3j. 

An edge-to-edge contact requires that the projection of the vector connecting the corre- 
sponding ends of two edges on the common perpendicular line of the two edges equals zero. 
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A nonoverlapping condition associated with an edge-to-edge contact requires that the pro- 
jection of the two edges on the common perpendicular hne of the two edges are completely 
separated. The 6 possible edge-to-edge contacts lead to 3 independent conditions, i.e., 

[vA - (vs - Ai + A2 + A3)] ■ lo > 0, (17) 

[vc - (vB - Ai + A3)] ■ li > 0, (18) 

[vc - (v,, + A3)] ■ I2 > 0, (19) 

where Iq = (va-V£;)x (vb-vd), h = (vb-v^) x (v^-vc) and I2 = (v^-v^) x (vc-v^). 
Since all of the quantities involved in Eqs. f ll7p -( ll9p are explicitly given, we do not provide 
their lengthy forms here. 

Furthermore, there are two additional nonoverlapping conditions given by 4 potential 
vertex-to-face contacts, i.e., 

A2-ni>-^| (20) 

where ni = (^, — ^^, — |) is the unit outward normal of the contacting face, 

(Ai - A2) ■ n2 > yi (21) 

and n2 = (— 3 , ~'q^-i \) is the unit outward normal of the contacting face. The edge-to- 
edge and vertex-to-face contacts are realized when the equality holds in the above conditions 
(iT5l) -(l2T]). Note we do not provide the lengthy forms of the above equations because the 
quantities involved have been given above. These nonoverlapping conditions constrain the 
possible values of the remaining three variables. 

Finally, the density of the six-parameter packings is given by 

where Vt = a/2/12 is the volume of a regular tetrahedron with unit edge length and Vp is 
the volume of the fundamental cell given by 
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Vf = |(a2/3i + ai/32)(-4y6 + Veai + Sv^az - 2^3 + 2^6/31 + /Jg) 

-l^(a3/3i + ai/33)(-2V3ai - lOas + Sv^ag + S/^a) (23) 

+^(-"3/32 + 02/33) (-6ai - 2v^«2 + V6a3 - 2x/6/33). 
The packing density relation (122|) in conjunction with (!23|) is a function of the 6 variables 
ttj and (3i {i = 1,2,3). We optimize the density for this problem using a sequential linear 
programming method subject to the nonoverlapping constraints given by Eqs. (1T5ll - (l2Tl) 
that specify all of the remaining possible contacts. In this way, we obtain the densest known 
packings with density (p = 0.856347. . .. This is a numerical verification of the optimality 
of the these packings among all four-particle basis packings. In the following section, we 
provide analytical constructions of the densest known packings and other dense packings 
from our general six-parameter family. 

B. Three-Parameter Family of Dense Tetrahedron Packings 

Both our simulations and local numerical analysis of Eqs. ( l22l) and ( l23l) suggest that 
realizing the face-to-face contacts associated with Eqs. ( ITSll and (fT6l) and the edge-to-edge 
contact associated with Eq. ( TT71) will lead to denser packings. Therefore, we let the equality 
hold in Eqs. f lTSj) . f fT6|) and f fT7|) . By solving these three additional equations, we can further 
eliminate three independent variables for the packing, i.e.. 



02 



V3/2, (24) 



/32 = -(-2v^ai + 2^2 + 372^3), (25) 

/33 = Y^i-^"^ + sVeai + 372^2 + 572^3 + IOVg/^i). (26) 

Therefore, the lattice vectors and the centroids of the dimers are functions of the three 
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parameters (ai,a3,/3i) only, i.e., 

_ 1 as 4 y/6 

^ f V^ ^ 

Ai = (-«!, — , -as), 

X t^ v^ , v^ 3V2 13^6 v^ 1 ;-„ , 

A2 = (A, - — + ^"1 ^«3, ^^^ ^tti - 2^3 - v6/:ii), 

. ,2 ^ 73 V3 3^2 21v^ ,76 1 

As = (--«!- Pi, -^ ^ai + ^pas, ^ + ^"1 - 2^3 + v6/3i). 

The density of the packing is given by 

_ 4VV _ V2 

where Vf is the vohime of the fundamental cell given by 



(27) 



(28) 



Vf = — (39V2 - 60V2al + 40\/3aias + 30V2al). (29) 

Note that although the packing structures depend on all three variables (ai,as,/3i), the 
density is only dependent on two variables (ai,as). The feasible values of the remaining 
variables are determined by Eqs. (TT7|) - (PT]) . In particular, we have 



3 


5V6 ^ 1 
12^^-20' 


3 

-ai- 


5^6 ^ 1 
" 12"^- 20' 



,7 VQ 7 VQ , n -rSl V633 ^6, 

max\ «! H as, as r < Pi < rmn\ — I — ai as, ai H as }■ 

^10 6 10 6 ^- ^4 2 4 ^' 4 2 4 ^^ 

(30) 

From Eq. fl25]l , we see that maximizing the density is equivalent to minimizing the volume 

of the fundamental cell Vp- To obtain our two-parameter family of packings reported in 

Ref . [36| , we assume that the minimum of Vp can be obtained by optimizing the two variables 

(ai,as) separately, i.e., we assume the following is true: 

min{y(^ai,a3)}VF = min{\/a^}[min{ya^}VF]. (31) 
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FIG. 2: (color online). The density (p surface of our family of tetrahedron packings as a function 
of the two parameters a and h (the parameters are related to the original variables via ai = a, 
as = — -v/|o, /?i = I — 2o — 6 for ai > and ai = a, as = — -v/|o, /3i = | + a — 6 for ai < 0). As 
explained in the text, the thick red lines (associated with a ~ ±0.016) show two sets of tetrahedron 
packings with distinct structures but with the same density. The two red points (associated with 
a = i]^) correspond to the densest two tetrahedron packings in this family. The packings found 



35l | are recovered from our two-parameter family (thick blue line, associated with 



by Kallus et al. 
a = 0). 

The inner optimization (over as with fixed ai) requires 

dVp/das = -{2V3ai + Sv^ag) = 





(32) 



which gives as = — ^/|ai. This further reduces the number of variables and leads to a 



two-parameter family of packings. Substituting this expression into Eq. ( 129|1 . we obtain 



VF = ^{39V2-80V2al). 

Substituting Eq. ( 133|) into Eq. (!28|) leads to the density expression 

100 



(33) 



(34) 



^ 117-240af 

where ai G (— jfo; ifo)- ^^ i^ important to note that for each ai 7^ 0, there are two sets 
of packings of tetrahedra, each with distinct structures but possessing the same density (as 
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shown in Fig. |2]by the thick red hnes). In addition, by substituting ai = a, as = — w|a, 
/3i = I — 2a — 6 for «! > and ai = a, a^ = — -i /|a, /3i = | + a — 6 for ai < 0, respectively 
into Eq. ( 1271) . the two sets of lattice vectors and dimer centroids can be obtained, i.e., for 



~i4o < ct < 0, we have 



1 a 4 3v^ [2 



Ai = (-a, — ^' Yg^)' 

7 /o^ Q -in 

A3 = (-- - 2a + 6, -- - V3a, -^ + -a - V66), 



(35) 



where < 6 < ^+3^°" : for < a < -^, we have 

dU ' 140 ' 



A a 4 3^2 [2 , 

^^ = ^5 + 3'-w!'~W!^Vr^' 

Ai = (-a, — ^' Y3«)' 

/3 ^ , ^3 ^ 3 10 ^,, 

As = (- - 2a - 6, -— + v/3a, -^ + 7|« + ^6&), 

A3 = (-- + a + 6,---v/3a,-^--^a-v/66), 



(36) 



where < 6 < 



3- 140a 
60 ■ 



The densest packings in this family are associated with a = —-^, h = Q and Ci = ^, 
6 = 0, possessing a density 0max = Ylfff = 0.855506.... In each set, there is a unique packing 
structure associated with (pmax (shown as the red points in Fig. [2]), instead of a spectrum of 
structures. Two different configurations of the densest packings of the two-parameter family 
within their corresponding rhombohedral fundamental cells are shown in Fig. [3]^a) and (b). 
The two four-particle configurations in the fundamental cells are only slightly different from 
one another. Specifically, the difference between the coordinates of the dimer centroids that 
are not at the origin (i.e., at Ci) is (f^,0, ^). Figure HJ^a) and (b) depict periodic rephcates 
of the densest tetrahedral packings corresponding to those shown in Fig. [3]^a) and (b) with 
8 fundamental cells (2 along each lattice vector). These packings are not chiral. The lattice 
vectors of packings associated with positive and negative a are related by an isometric map 



37|. 
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(a) 



(b) 





(c) 



(d) 



FIG. 3: (color online). In the top panel [(a) and (b)], we show the two different configurations of 
the densest packings of four tetrahedra (two dimers) with (j) = ^^^ = 0.855506 . . . within their 
corresponding rhombohedral fundamental cells of the two-parameter family, (a) Here a = — ^^ 
and 6 = 0. (b) Here a = ^ and 6 = 0. The fundamental cells are colored yellow (light gray) and 
their boundaries are colored red (dark gray). The two packings and their associated fundamental 
cells shown in (a) and (b) are only slightly different from one another. In the bottom panel [(c) 
and (d)], we show the two different configurations of the densest packings of four tetrahedra (two 



dimers) with cf) = ^g^j = 0.856347 . . . within their corresponding rhombohedral fundamental cells 
of the three-parameter family, (c) Here ai = j^, as = jgoV I ^^^^ Pi = xgo' (d) Here ai = — jggj 



03 = "iloV 3 ^'^d /3i = §5- The fundamental cells are colored green (medium gray) and their 
boundaries are colored red (dark gray). The two packings and their associated fundamental cells 
shown in (c) and (d) are only slightly different from one another. 
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(a) 



(b) 





(c) 



(d) 



FIG. 4: (color online). This figure shows periodic replicates the densest tetrahedron packings of 
the two-parameter family (with (j) = j^gjg = 0.855506. . .) and the three-parameter family ( with 
(f) = |§^ = 0.856347...) corresponding to the top and bottom panels in Fig. El respectively, 
with 8 fundamental cells (2 along each lattice vector). The tetrahedra within the fundamental 
cells are shown in blue and green (both appearing as medium gray in the print version) for the 
two-parameter [(a) and (b)] and three-parameter [(c) and (d)] family packings, respectively. Note 
that in (a), the dimer with centroid at Ai (blue, or medium gray in the print version) is slightly 
shifted to the right with respect to the dimer at the origin O (yellow, or light gray in the print 
version); and in (b), the dimer with the centroid at Ai (blue, or medium gray in the print version) 
is slightly shifted to the left with respect to the dimer at the origin O (yellow, or light gray in 
the print version). The three-parameter packing structures are very similar to that of the two- 
parameter packings. Observe that in the perspective shown the dimers appear as if they were 
single tetrahedra. 
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At a = 0, there is only one set of packings with the same (j) = -^ = 0.85470..., whose 
structures are dependent on h (shown as the blue line in Fig. [2]). These packings reduce 
exactly to those discovered by Kallus et al., which possess relatively high symmetry. In 
particular, a = allows the centroids of the dimers, which are related to each other by 
an integer multiple of Ai, to be perfectly aligned on two of the mirror image planes of 
the dimers simultaneously, which leads to additional two-fold rotational symmetry of the 
packing. This additional rotational symmetry, together with the point inversion symmetry, 
leads to uniform packings with respect to each tetrahedron (not just each dimer), i.e., the 
symmetry operation acts to take each tetrahedron into another. 

If we do not make the assumption used above (i.e., optimizing the two variables ai and 
as separately), the boundary of the domain of Vp needs to be searched completely to find 
the minimal value. For the boundary specified by 

3 5v^ 1 

-Q!i a^ < — , (37) 

2 12 20 ^ ' 

we have ai = ^(3 + 25-\/6q;3). Substituting this expression into Eq. (I29l) gives 

Vf = —— (219^2 - sVSaa + 20(}V2al) (38) 

the minimal value of which is Vp = „iff'^ with a^ = yIo \/ f- -^°^ ^^^ boundary specified by 

3 5v^ 1 

-Oil as > , (39) 

2 12 ^ - 20' ^ ^ 

we have ai = ^(—3 + 25\/6as). Substituting this expression into Eq. fl29|) gives 

Vf = -^(219^2 + sVSas + 200V2al) (40) 

the minimal value of which is Vp = gnoo/o "with a^ = ■"yIo V §■ 

Thus, we find a maximal density (pmax = j^, which was also reported by Chen et ai 



37(1, associated with ai = tIqj'^s = tIoa/s ^^^ '^i ~ "~T60''^3 ~ "leoA/s- '^^^ packing 



160' ■^ 160 V 3 ^ 160' ■^ 160 Y 3' 

3^^. + via, = iii 



structures can be obtained by substituting /3i = | — |ai + ^as = ^ ior ai > and 
/3i = I + |ai — ^as = II for ai < into Eq. (p7|) . These three-parameter packings are 
very close to the two-parameter packings in structure and possess the same symmetry [see 
Fig.[3](c) and (d), and Fig. SJ^c) and (d)]. 
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III. TOWARDS UPPER BOUNDS ON THE MAXIMAL DENSITY 

The problem of determining upper bounds on the maximal density of packings of non- 
spherical particles is highly nontrivial, and yet such estimates would be indispensable in as- 
sessing the packing efficiency of a candidate dense packing, especially if tight upper bounds 
could be constructed. It has recently been shown that (pmax of a packing of congruent 
nonspherical particles of volume vp in M? is bounded from above according to 

Vp 71 



^max < min 



pi 



(41) 



where vs is the volume of the largest sphere that can be inscribed in the nonspherical particle 



and 7r/vl8 is the maximal sphere-packing density |28|, |29|. The upper bound fj4T]) will be 
relatively tight for packings of nonspherical particles provided that the asphericity 7 (equal 
to the ratio of the circumradius to the inradius) of the particle is not large. However, for 
tetrahedra, the asphericity is too large for the upper bound fHT]) to yield a result that is less 
than unity. 

One possible approach to obtaining nontrivial upper bounds is to attempt to generalize 



the idea that Rogers used to prove upper bounds on 4>max for sphere packings [38[. The 
key concept is to consider a locally dense cluster of 4 contacting spheres in a tetrahedral 
arrangement and then prove that the fraction of space covered by the spheres within the 
tetrahedron joining the sphere centers is an upper bound on (pmax- This can be done because 
one can triangulate any sphere packing to decompose it into generally irregular tetrahedra 
with vertices at sphere centers. The fact that the regular tetrahedron has the best density 
for any tetrahedron, then yields an upper bound for the density of any sphere packing. In 
the case of the non-tiling Platonic and Archimedean solids, a natural choice for the enclosing 
region associated with the cluster is its convex hull. 

For tetrahedra, we must identify the densest local cluster with density that exceeds or 
equals (pmax while taking into account their shape symmetry. A trivial choice is a dimer 
because the fraction of space covered by the dimer within its convex hull is unity. Two 
nontrivial choices are the 5-particle "wagon-wheel" cluster as shown in Fig. EJ^a) and a 
more "isotropic" icosahedron-like cluster (a cluster consisting of 20 tetrahedra that share a 
common vertex) as shown in Fig. [5t^b). The "wagon- wheel" and icosahedron-like cluster can 
be considered to be the most anisotropic and isotropic densest local clusters, respectively. 
The former is essentially a "fiat" (quasi-two-dimensional) structure. Since neither of the 
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(a) 



(b) 



FIG. 5: (color online), (a) The convex hull of five regular tetrahedra in a "wagon wheel" arrange- 
ment. The convex hull can be decomposed into five regular tetrahedra (shown in blue, or medium 
gray in the print version) and five thin irregular tetrahedra (shown in red, or dark gray in the print 
version). (b)The convex hull of 20 regular tetrahedra in an "icosahedral" arrangement. The convex 
hull can be decomposed into 20 regular tetrahedra (shown in blue, or medium gray in the print 
version), 12 pyramids with pentagonal bases (shown in yellow, or light gray in the print version) 
and 30 pyramid with rectangular bases (shown in red, or dark gray in the print version). 



two local clusters can tile space, their local packing densities, which are computed below, 
provide limiting estimates on upper bounds for 4>max- 

The convex hull of the "wagon-wheel" cluster composed of five tetrahedra sharing a 
common edge (a "wagon wheel" cluster) can be decomposed into five regular tetrahedra 
and five thin irregular tetrahedra. We assume the gaps between the regular tetrahedra are 
equal, i.e, the thin irregular tetrahedra are congruent. Note the regular tetrahedron shares 
two faces with its neighboring irregular tetrahedra. Thus, the volume ratio is equal to the 
ratio of the corresponding heights of the regular and irregular tetrahedron, i.e., 

Vt 1 



7 



(42) 



Vn ^^(3cos2ff-l)' 

where Vr and Vr, is the volume of the regular and irregular tetrahedron, respectively. Thus, 
the density of this local packing, defined as the fraction space covered by the regular tetra- 
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hedra within the convex hull is given by 

Vt 



>w 



Vt + Vt, 3^cOs2f^ + (l-^) 



0.974857.... (43) 



Because this cluster is highly anisotropic and is an effectively "ffat" object, it is reasonable 
to assume that it is not the least densest local cluster and therefore its local density is likely 
to be gross overestimate of the best upper bound on (pmax- 

It seems reasonable to suppose that the least densest local cluster with the highest shape 
symmetry is one consisting of 20 tetrahedra sharing a common vertex, i.e., an icosahedron- 
like cluster. The convex hull of such a cluster can be decomposed into 20 regular tetrahedra, 
12 pyramids with pentagonal bases and 30 pyramid with rectangular bases [see Fig. EJ^b)]. 
We assume the gaps between tetrahedra are equal and thus the two types of pyramids are 
congruent. The volume of the regular tetrahedron is Vt = ^ and the volume Vp of the 
pyramid with pentagonal base is given by 



Vp = ltan(— )L%/l ^, (44) 

12 UO^ Y 4cos2ff' ^ ' 



and the volume Vr of the pyramid with rectangular base is given by 




Vr = ^LvO^L^, (45) 



where 

and r = (1 + -\/5)/2 is the golden ratio. Thus, the local packing density is given by 

(t)i = '^^ = 0.880755 . . .. (47) 

^ 20l^r + 12Vp + 30Fr ^ ^ 

Note in the above calculations, we have assumed that the gaps between the tetrahedra are 
equal, which is sufficient to provide an estimate of the fraction of space covered by the cluster 
within its convex hull. It is noteworthy that the icosahedron-like cluster of 20 tetrahedra 
is at least a local extremum, i.e., adding or removing a tetrahedron from the cluster will 
decrease the local packing density. Moreover, adding another "shell" of tetrahedra around 
the icosahedron-like cluster will also reduce the local packing density. We note that the idea 
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of Rogers to prove an upper bound for sphere packings cannot be used here because there is 
no analogous decomposition of space into irregular convex hulls shown in Fig. |5t^a) or (b). A 
completely new idea is needed to prove that the aforementioned estimates are bounds. If one 
could prove that these clusters are indeed locally denser than the globally densest packings, 
then the aforementioned estimates provide upper bounds on (pmax, but they cannot be sharp, 
i.e., they are not achievable by any packings. Without loss of generality, it is reasonable to 
expect that the globally densest packings of tetrahedra contain a particular proportion of 
both the isotropic and anisotropic local clusters. Therefore, these limiting estimates provide 
a plausible range of upper bounds, i.e., (p^ax ^ [0.880755, 0.974857]. It is noteworthy that 
the density = 0.856347. . . of the densest known packings of tetrahedra is relatively close 
to this putative lower-limit upper bound density estimate of 0.880755 . . .. 

IV. CONCLUSIONS AND DISCUSSION 

For all of the small periodic packings that we investigated (including 2 to 32 particles 
per fundamental cell) using our numerical ASC scheme, the densest packings that emerged 
had a 4-particle basis. The numerical optimization that we performed of our six-parameter 
family of tetrahedron packings verified the optimality of the densest 4-particle packing with 
density (p = Ip^ — 0.856347.... Are these dimer packings optimal among all packings? 
This is a difficult question to answer in any definitive way. We first note that the two con- 
tacting dimers in the fundamental cell possess center inversion symmetry and thus they can 
be viewed as a centrally symmetric compound object. Hence, our two-parameter family of 
packings for the two-dimers (and the densest known packings) are Bravais lattice packings of 
such compound objects. Thus, it is not very surprising that dense Bravais lattice packings of 
this centrally symmetric compound object have a fairly high density based on the arguments 



leading to Conjectures 1 and 2 [28|, |29|, i.e., there are a large number of face-to-face contacts 
which are made possible due to the central symmetry of the compound object and bring the 
centroids of the objects closer to each other. However, Conjecture 3 states that central sym- 
metry alone does not guarantee optimality of the densest Bravais lattice packings among 
all packings. That is, there are additional constraints on the geometry and shape of the 
objects, e.g., they might also need to possess three equivalent principal axes, as do the cen- 
trally symmetric Platonic and Archimedean solids. When the latter condition does not hold. 
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there could exist non-Bravais lattice packings that are denser than the optimal lattice pack- 



ings (e.g., ellipsoids |2l|). The 4-tetrahedron compound object is very anisotropic and does 



not possess equivalent principal axes. Moreover, it is a concave object rather than convex. 



which makes it even more difficult to make a definitive cone 



usion about the optimality of the 



densest dimer packings. In addition, Haji-Akbari et al. J3J] have found periodic packings of 
tetrahedra consisting of 82 particles per fundamental cell with a slightly lower density (i.e., 
(f) ~ 0.8503) than that of the densest dimer packings. Their packings possess much larger 
fundamental cells with distinctly different and more complicated particle arrangements than 
the dimer packings. This suggests that there is a large degeneracy of different non-Bravais 
lattice packing structures with high densities with various levels of complexity and perhaps 
even higher densities than the optimal 4-particle packings. Such denser packings could be 
discovered by carrying out exhaustive searches to determine the globally maximal densities 
of packings with successively larger numbers of particles per fundamental cell. 

Previous numerical studies have indicated that dense packings may have a large number 
of particles per fundamental cell arranged in a complex fashion, e.g., the "disordered wagon- 



wheels" packing with = 0.822637 . . . (29[ and the "ring stacks" packing with = 0.8503 . . . 
34j |. However, it is now clear that such packings are in fact only locally optimal solutions and 
hence the numerical techniques used to obtain them are incapable of extricating themselves 
from these "trapped" regions of configuration space to find denser and more ordered struc- 
tures due to the intrinsic geometrical frustration of the tetrahedron mentioned earlier. This 

n 

may also call into question claims made by Haji-Akbari et al. |34] that their packings, which 
are characterized by an effective "quasicrystal-like" plane, are true thermodynamic equilib- 
rium phases of tetrahedra at high densities. Instead, our high-density constructions suggest 
that uniform periodic packings with a 4-particle basis (or even some yet unknown denser pe- 
riodic packing) and their unjammmed, lower-density counterparts could be the stable phases 
at high densities. If the latter is correct and the putative "quasicrystal-like" phase truly ex- 
ists at intermediate densities, then it is hard to imagine how such complex quasicrystal-like 
structures of tetrahedra under quasi-static compression (densification) could rearrange them- 
selves at higher densities to spontaneously form a more ordered periodic arrangement with 
higher symmetry. However, it is difficult to draw any such definitive conclusions without 
further study. 

Although there could still be tetrahedron packings denser than our constructions, it ap- 
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pears that all of the evidence thus far points 



tot 



cannot possess very high symmetry 28 



-30 



33 



le fact that the densest tetrahedron packings 



34j due to the lack of central symmetry of a 



tetrahedron and because tetrahedra cannot tile space [29]. Indeed, the dense 4-particle-basis 
sackings found in the present paper improved upon the best dimer packings of Kallus et al. 



35| by relaxing the two-fold rotational symmetry constraints they imposed. 
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